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Teleportation is the most widely discussed application of the basic principles of quantum mechan¬ 
ics. Fundamentally, this process describes the transmission of information, which involves transport 
of neither matter nor energy. The implicit assumption, however, is that this scheme is of inherently 
nonlocal nature, and therefore exclusive to quantum systems. Here, we show that the concept can 
be readily generalized beyond the quantum realm. We present an optical implementation of the 
teleportation protocol solely based on classical entanglement between spatial and modal degrees of 
freedom, entirely independent of nonlocality. Our findings could enable novel methods for distribut¬ 
ing information between different transmission channels and may provide the means to leverage 
the advantages of both quantum and classical systems to create a robust hybrid communication 
infrastructure. 


INTRODUCTION 

Quantum teleportation allows a sender to share in¬ 
formation with a receiver by encoding it into a quan¬ 
tum state that both parties can access. In general, the 
counter-intuitive aspects of this type of protocols stem 
from the fact that, in contrast to all conventional forms of 
communication, the state itself is not being transmitted 
between the two [1] . Following the first landmark experi¬ 
ments employing photonic quantum states as information 
carriers mu, teleportation has been implemented in a 
wide variety of different physical systems ISHH]- Telepor¬ 
tation is therefore considered to be a fundamental build¬ 
ing block for quantum communication and computation 
architectures. At the core of its ability to distribute in¬ 
formation through potentially vast and far-reaching net¬ 
works is entanglement, which allows the nonlocal inter¬ 
action that is necessary for a transfer of information to 
occur [9]. 

For many years, entanglement was viewed as the epit¬ 
ome of a quantum effect, seemingly lacking any analogue 
in classical systems. Only very recently it became clear 
that the underlying algebraic structure giving rise to en¬ 
tanglement can indeed be reproduced in the classical do¬ 
main [MflE]. Very recently, moreover, state transfer pro¬ 
tocols based on the classical entanglement concept have 
also been proposed m- As it turns out, the key distinc¬ 
tion between classical and quantum systems is not their 
capability to support entanglement, but rather the (non- 
) locality of the same m- As subsequent works have 
demonstrated [THHH, these notions may help to infuse 
classical physics with the potentialities offered by quan¬ 
tum entanglement to drive novel and interesting appli¬ 
cations. In this work, we extend the concept of telepor¬ 


tation beyond its well-known quantum context. To this 
end, we propose and experimentally demonstrate a pho¬ 
tonic setting in which classical entanglement can be har¬ 
nessed to teleport information between the spatial and 
modal degrees of freedom of a purely classical light field. 
While nonlocality remains the differentiating feature be¬ 
tween quantum and classical configurations, our results 
illustrate that teleportation-based protocols can bridge 
the gap between these two realms. 

THEORETICAL BACKGROUND 

On a conceptual level, teleportation involves at least 
three degrees of freedom, which henceforth will be la¬ 
belled as A, B and C. Whereas an entangled Bell state is 
established between A and B, the information one desires 
to transmit is contained in C. The actual teleportation 
process is facilitated by the application of a controlled- 
NOT (C-NOT) operation with respect to A and C, fol¬ 
lowed by a Hadamard operation and the subsequent pro¬ 
jection onto a Bell state [7]. For the case of quantum 
teleportation, these degrees of freedom are represented 
by qubits, while in the classical case they are called cebits 
[16] . By adopting Diracs notation, the input state can be 
written in the following form (see Supplementary Mate¬ 
rial) 

1^) = ^ [(«io>c+ me) ® , (1) 

where the cebit C resides in a two-dimensional Hilbert 
space. The complex coefficients a and /3 carrying the 
information can be normalized such that \a\‘^ + = 
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1. Here, \^~^)ab represents a classically entangled Bell 
state. The aforementioned combination of C-NOT and 
Hadamard gates then yields the following output state: 
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The entanglement, initially shared between the cebits 
A and B, has been transferred to C and A. To realize the 
teleportation of information between the cebits C and 
B, I Tout) is projected onto one of the four basis states 
spanning the joint Hilbert space of the two entangled 
cebits A and C. A readout of the information contained 
in the cebit B will therefore give us access to the informa¬ 
tion initially encoded in the cebit C. It is worth noting 
that Eq. <© is completely independent from the charac¬ 
ter of the degrees of freedom involved. Rather, it is the 
choice of entanglement that contextualizes the associated 
process in the classical, local, world, or in the nonlocal 
realm of quantum mechanics. This observation leads us 
to interpret Eq. <§ as the teleportation equation. It di¬ 
rectly follows that the gate sequence of the teleportation 
scheme is universally valid, irrespectively on the nature 
of the entanglement employed. 


EXPERIMENT 

To implement such a protocol, we employ the exper¬ 
imental setup shown in Eig. where the information 
can be teleported from the path (C) to the spatial (B) 
Hilbert space. Light from a 594 nm He-Ne laser is sent 
to a custom- designed wave plate (RPWP) to synthetize 
a single classically entangled beam with radial polariza¬ 
tion, corresponding to \^~^)ab [ 12 ] • Having prepared the 
cebits A and B for our classical implementation of tele¬ 
portation, the third cebit C can be readily established 
between two identical copies of such classically entan¬ 
gled beam, propagating along different paths. One path 
represents |0)c, and the other \l)c^ respectively [14]. Ex¬ 
perimentally, this is realized by means of a 50/50 beam 
splitter (BSl), and two neutral-density filters {F^ and 
F^) placed at each output channel of BSl encode the 
information in the cebit C. 

The output state |Tout) described by Eq. (© therefore 
corresponds to the two output channels of BS2. Note 
that each of its four unique mode structures can be indi¬ 
vidually extracted with a simple measurement operation. 
Consider, for example, the first term in Eq. § , which 
can be written, in terms of optical modes, as follows (see 
Supplementary Material): 



FIG. 1. (a) Generation: A continuous-wave He-Ne laser is 
prepared in a classically entangled state with radial polar¬ 
ization by a rotating polarization wave plate (RPWP). The 
initial state |T) is then obtained using a beam splitter (BSl). 
Two density filters (Fh and Fp) are used to encode the infor¬ 
mation in the cebit G. Preparation: The initial state |T) is 
then sent through a G-NOT gate, realized using two Sagnac 
interferometers with a polarizing beam splitter. The lower in¬ 
terferometer, corresponding to ll)^ also contains a half-wave 
plate (HWP) to rotate the polarization. The two parts of 
the beam are then recombined through a second beam split¬ 
ter (BS2), which implements a Hadamard operation for the 
cebit G. Measurement: The correct output state is selected by 
choosing the x-polarization (PBS3) of the lower output chan¬ 
nel of BS2. The reflected beam from a third beam splitter 
(BS3) is sent to a GGD camera (Gam2) for direct acquisition 
of the intensity profile and angle measurement. The transmit¬ 
ted beam from BS3 is instead sent to the modal decomposition 
stage, consisting of a computer generated hologram (GGH) 
a lens and a second GGD camera (Garni), (b) Exemplary 
scanning electron microscope image of a nanograting-based 
polarization rotating wave plate such as the RPWP used to 
generate the classically entangled beam with radial polariza¬ 
tion. 


|00)cA G) (q^|0)b + P\^)b) 

= F p'lpoiix.y)]. (3) 

Therein, t/iq and t/qi are the first-order Hermite- 
Gaussian solutions of the paraxial equation [22], and 
x,y denote the horizontal and vertical polarization eigen¬ 
states, respectively. The path cebit C has been imple¬ 
mented by labelling the upper and lower beam with the 
superscript u and /, respectively. Then,|0)c' corresponds 
to the beam being in the upper path, while \l)c corre¬ 
sponds to the beam being being in the lower path [T4] . 
A closer inspection to the above equation reveals that in 
order to project the output state |Tout) onto the state 
it is sufficient to take the x-polarization eigen¬ 
state of the beam in the upper path, i.e., |0)c- This 
guarantees access to the information encoded in the cebit 
B. In our experimental setup, this projection operation 
is realized by placing a polarizer (PBS3) at the output 
of BS2 corresponding to |0)c- In terms of optical beams. 
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FIG. 2. Estimation of the ratio ^ j ol initially encoded in the 
cebit C as retrieved from cebit B by measuring the angle 
formed by the intensity distribution acquired from the CCD 
camera Cam2 with the vertical axis. As can be seen, the re¬ 
trieved angles (black points) are in good agreement with the 
theoretical prediction (red solid line), namely cot^ = |/3a|. 
The insets show the acquired intensity distribution (single re¬ 
alization) corresponding to^ = 47°, 55°, 62° and 76°, respec¬ 
tively. 


the projected state cA(00|^out) assumes the following 
form: 

M{x,y) = - yo) + - yo). (4) 

The beam M{x,y) is then split into two parts by a 
beam splitter (BS3). While the reflected beam is directly 
acquired by a CCD camera (Cam2), the transmitted 
beam is instead sent through a computer-generated holo¬ 
gram (CGH) that decomposes it in its different Hermite- 
Gaussian constituents [23]. Subsequent imaging onto a 
second CCD camera (Caml) allows for the extraction of 
both the amplitude and the phase of the modal coeffi¬ 
cients a and /3. 

The information encoded in the spatial cebit B can 


be then retrieved in two different ways. If the phase 
difference between a and /3 is negligible, the intensity 
distribution acquired by Cam2 is a Her mite-Gaussian 
mode '001 rotated by an angle 0 with respect to the y- 
axis. The value of the modal coefficients can be then 
directly obtained from the acquired image by measuring 
the rotation angle 0 with a digital algorithm (see Meth¬ 
ods). Exemplary results of this measurement, together 
with the acquired intensity distributions corresponding 
to 0 = 47°, 55°, 62° and 76°, are shown in Fig. and 
follow the theoretical prediction \Pa\ = cotO. If, on the 
other hand, the phase difference between a and [3 is not 
negligible, and a full modal decomposition is required to 
extract the information, as illustrated in Fig. 

Clearly, the choice of the spatial degree of freedom as 
the target cebit for teleportation allows for an easy and 
reliable measurement and a classical communication net¬ 
work, thus opening the possibility for new ways to reroute 
information between initially disjoint Hilbert spaces. 

Generation of the classically entangled state 

To generate a classically entangled beam with radial 
polarization, we employed a rotating polarization wave 
plate (RPWP), like the one shown in Fig. Sb). This 
special optical component was custom-fabricated using 
a femtosecond laser writing technique and it consists 
of a A2 plate with locally varying anisotropy. The lo¬ 
cal anisotropy is realized by inscribing periodic sub¬ 
wavelength nano-gratings in the bulk of a fused silica 
substrate [2S Eg. By tuning the polarization of the 
inscribing laser, the orientation of the nanograting can 
be set as desired. This offers a great degree of control 
over the polarization pattern imprinted on the impinging 
beam on the micron scale. 


Estimation of the modal coefficients via angle 
measurements 

When the phase difference between a and /3 is negligi¬ 
ble, the value of these modal coefficients can be retrieved 
by directly acquiring the modal function M{x^y). Ac¬ 
cording to Eq. in fact, the intensity distribution 

acquired by Cam2 consists of a Hermite-Gaussian mode 
-001 rotated by an angle 0 with respect to the vertical 
axis. To measure such an angle, we first evaluate numer¬ 
ically, from the retrieved image, the line L joining the 
center of mass of each lobe. Then, we compute the slope 
of the line perpendicular to L, to directly obtain the an¬ 
gle formed by the intensity distribution |M(x,^)p with 
the vertical axis. The values of a and p can be therefore 
retrieved by using the relation cot?^ = \Pa\ (see Sup¬ 
plementary Material). In order to provide an accurate 
measurement of the angle the results presented in Fig. 
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b) 

A0 = 0.1822 ±0.0071 A0 = 1.4909 ± 0.0539 A(^ = 3.829 ± 0.0596 



FIG. 3. (a) Plot of the ratio \P/a\ as retrieved from cebit B 
as a function of the ratio |/3/(a| initially encoded into cebit C. 
(b) Comparison of the retrieved intensity distributions (upper 
row) and the theoretical predictions (lower row) correspond¬ 
ing to different values of the phase difference A0 encoded into 
a and /3 as retrieved from cebit B. To encode the phase in¬ 
formation into the cebit C, the filter Fa has been tilted by 
2.5°, 5° and 10°, respectively, while the filter Fp has been left 
unchanged. The correspondent retrieved phases are displayed 
above the acquired intensity distributions, together with the 
corresponding error intervals. 


[^have been averaged over fifteen acquired intensity dis¬ 
tributions corresponding to the same set of parameters 
a and /3. The error bars in Fig. moreover, correspond 
to the standard deviation calculated with such a set of 
experimental data. 


to retrieve the information encoded in the spatial cebit 
B. With this technique, the modal function M{x,y) de¬ 
scribed by Eq. @ is decomposed into the Hermite- 
Gaussian modes ipiQ and '0oi- To do that, we employ the 
correlation filter method (CFM) technique (see Supple¬ 
mentary Material) [23], which basically implements the 
scalar product between the modal function described by 
Eq. and the Hermite-Gaussian basis function 
i.e., 

= f f 'ipl^{x,y)M{x,y)dxdy. (5) 

J J —oo 

Therefore, since the Hermite-Gaussian functions 2 pio 
and 2 poi are the complete basis set that spans the two 
dimensional Hilbert space associated to cebit B, this 
method allows to extract the information encoded in this 
cebit by simply calculating ('0oi|M) = a and ('0io|M) = 

respectively. This is realized in our experimental setup 
by using a computer-generated hologram (CGH) with ap¬ 
propriately designed transmission functions [23| . The in¬ 
tegration appearing in Eq. (§ has been then carried 
out using a 2/ system, which performs the Eourier trans¬ 
form of the incoming signal. The intensity distribution 
in the focal plane therefore provides access to the ampli¬ 
tude of the coefficient a and With this method, it is 
also possible to retrieve the relative phase difference A0 
between a and P by means of a three-way interference be¬ 
tween different replicas of the two mode functions with 
known phase difference. In addition, in order to provide 
real-time access to both a and [3 simultaneously, we em¬ 
ployed an angular multiplexing technique similar to the 
one described in Ref. [23] . 

Although the results presented in Eig. are in very 
good agreement with the expected ones, it has to be 
noted that the measurement error increases with the 
measured ratio [3/a. To obtain a good correlation sig¬ 
nal, in fact, the scalar product (§ needs to be evaluated 
at a single point (e.g., along the optical axis of the sys¬ 
tem) , where the signal originating from the two Hermite- 
Gaussian modes are truly independent of one another. 
Eluctuations of the position of this point due to small 
fluctuations of the position where the two beams recom¬ 
bine (BS2) are the main source of measurement errors 
and can be minimized by ensuring a sufficient contrast 
in signal strengths between the modes. 


CONCLUSIONS 


Estimation of the modal coefficients by modal 
decomposition 

When the phase difference between a and [3 is not neg¬ 
ligible, a full modal decomposition technique is needed 


With this work we have extended the concept of tele¬ 
portation to the classical realm by implementing a pho¬ 
tonic setting in which classical entanglement allows to 
teleport information between path and spatial degrees of 
freedom of a purely classical light field. In doing so, we 
have shown that teleportation is a general concept, that 
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transcends the distinction between classical or quantum 
systems, and that non-locality ultimately differentiates 
between these two realms. Along these lines, our results 
may pave the way towards the realization of a hybrid 
classical-quantum communication infrastructure. 
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